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ABSTRACT 

We combine two complete, X-ray flux-limited surveys, the ROSAT Bright Survey 
(RBS) and the ROSAT North Ecliptic Pole (NEP) survey, to measure the space 
density (p) and X-ray luminosity function ($) of non-magnetic CVs. The combined 
survey has a flux limit of Fx > 1-1 X lCP 12 ergcm _2 s _1 over most of its solid angle 
of just over 2n, but is as deep as ~ 10 _14 ergcm _2 s _1 over a small area. The CV 
sample that we construct from these two surveys contains 20 non-magnetic systems. 
We carefully include all sources of statistical error in calculating p and by using 
Monte Carlo simulations; the most important uncertainty proves to be the often large 
errors in distances estimates. If we assume that the 20 CVs in the combined RBS and 
NEP survey sample are representative of the intrinsic population, the space density 
of non-magnetic CVs is x 10~ 6 pc -3 . We discuss the difficulty in measuring $ in 
some detail — in order to account for biases in the measurement, we have to adopt a 
functional form for $. Assuming that the X-ray luminosity function of non-magnetic 
CVs is a truncated power law, we constrain the power law index to —0.80 ± 0.05. It 
seems likely that the two surveys have failed to detect a large, faint population of 
short-period CVs, and that the true space density may well be a factor of 2 or 3 larger 
than what we have measured; this is possible, even if we only allow for undetected 
CVs to have X-ray luminosities in the narrow range 28.7 < log(Lx/ergs _1 ) < 29.7. 
However, p as high as 2 x 10~ 4 pc -3 would require that the majority of CVs has X-ray 
luminosities below Lx = 4 x 10 28 ergs _1 in the 0.5-2.0 keV band. 

Key words: binaries - stars: dwarf novae - novae, cataclysmic variables - X-rays: 
binaries - methods: observational, statistical. 



1 INTRODUCTION 

Cataclysmic variables (CVs) are interacting binary stars 
consisting of white dwarfs accreting from low-mass, Roche 
lobe filling companions fsee lWarnerlll995l for a review). An- 
gular momentum is continuously lost from the binary orbit, 
driving mass transfer, as well as evolution in orbital period 

(Port). 

There are still many uncertainties in the theoret- 
ical descr i ption of CV formation and evolution (e.g . 
Kolbll2002l: llvanova fc Taaml |2003l ; iNelemans fc Toul 120051 ; 



HSB^ l ^ serious disc r^- 
cies between the predictions of theory and the proper- 
ties of the observed CV population (e.g. Patterson! 1 19981 : 
lAungweroiwit et al.ll2006l : iPretorius et al. I l2007al : Pretorius 
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& Knigge 2008a,b; Gansi cke et al 1 l2009l : lKnigge~ et al ] |201ll l. 
In order to constrain evolution models, more and better ob- 
servational constraints on the properties of the Galactic CV 
population are needed. A fundamental parameter predicted 
by CV evolution theory, that is expected to be more easily 
measured than most properties of the intrinsic CV popu- 
lation, is the space density, p. The luminosity function ($) 
is a closely related property; although more challenging to 
constrain than p, $ contains information on the mass trans- 
fer rate (M) distribution of CVs and is therefore potentially 
much more valuable than the space density alone. 

Some theoretically predicted v alues of the C V space 
density ar e (0.5 - 2) x 10 -4 pc -3 (|de Kooll fl992l). 1.8 x 
10" 4 pc -3 (|Kolbll 19931 ). 2 x 10" 5 pc" 3 (|Politandll99r3 ). Ob- 
servational estimates are typically lower, but have a large 
range; values from ^ 5 x 10~ 7 pc~ 3 to p ~ 10~ 4 pc~ 3 have 
been reported (e.g. | Ritter fc Burkertl[l986l ; iRitter fc Ozkanl 
Il986l : iHertz et al.) [l990l : IShara et al.lll993l : |Pattersonlll998l : 
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Warner 2001; Schreiber & Gansicke 


120031; Cieslinski et al. 


2003 


; lArauio-Betancor et al.l 2005 


; Aungweroiwit et al. 


2006 


; Rogcl ct al. 2003). Fewer estimates of the CV luminos- 



Bvckling et alj |2010| ) and hard X-ray i Revnivtsev et all 
20081) luminos ity functions have been published (see also 



Ak et al1l2008 | for near-IR and opt ical luminosity functions) . 



The work of iBvckling et al.l l|2010l ) is based on the sample of 
non-magnetic CVs with parallax distances less than 200 pc; 
although this sample is certainly not complete, it gives a 
firm low er limit on the lumino sity f unction. The CV samples 
used by ISazonov etafl l|2006h and IRevnivtsev et all (|2008h 
consist mostly of magnetic CVs. 

Uncertainty in p (and $) measurements is in part 
caused by statistical errors, arising from uncertain distances 
and small number statistics. However, the dominant source 
of uncertainty, as well as the cause of inconsistencies between 
some estimates, is most likely systematic errors caused by 
selection effects. The selection effects acting on CV samples 
are most easily accounted for in samples with simple, well- 
defined selection criteria. 

Whereas optical CV samples always include selection 
criteria based on, e.g., colour, variability, or emission lines, 
there are X-ray selected CV samples that are purely flux- 
limited. The completeness of these samples is easier to model 
than that of any existing optically selected sample. Further- 
more, there is a well-known empirical correlation between 
the ratio of optical to X-ray flux (F op t/Fx) and optical lu- 
minosity, implying that an X-ray flux limit does not intro- 
duce as strong a bias against intrinsi cally faint, short-period 
CVs a s an optical flux limit (e.g. Patterson &; Raymond! 
19854 l yan Teeseling fc Verbuntl 1 1994 Ivan Teeseling et al.1 
19961 ; iRichmanll 199rJ ) ~ 

The ROSAT Bright Survey (RBS) consists of all bright 
(count rate > 0.2 s _1 ), high Galactic latitude (|6| > 
30°) source s in t he R OSAT All-Sky Su rvey (RASS; see 
IVoges et al.1 Il999l and IVoges et al.1 |2000D . Optical coun- 
terparts have been i dentified for all sources in the RBS 
i|Schwope et all l2000h . and it inclu des 46 CVs (of whic h 
1 1 were previously unknown; see ISchwope et al. I l2002t ). 
ISchwope etail l|2002l ) already used this survey to estimate 
p ~ 3 x 10 -5 pc -3 for non-magnetic CVs; they note that 
this measurement is dominated by two systems with very 
small distance estimates, and would be an over-estimate 
if these two distances where sign ificantly under-estimate d 
(as has turned out to be the case; iThorstensen et al.i r2006). 
Omitting those two sources, ISchwope et al. (2002) obtain 
p= 1.5 x 10~ 6 pc- 3 . 

In a previous paper (jPretorius et al.1 l2007bT l. we have 
considered non-magnetic CVs from the ROSAT North 
Ecliptic Pole (NEP) survey; this survey is 2 orders of magni- 
tude deeper than the RBS, but covers a much smaller area 
(~81 sq.deg., down to roughly 10~ 14 ergcm~ 2 s~ 1 ; see e.g. 
iGioia et al.ll2003l ; lHenrv et aljfeood ). and also has complete 
optical follow-up. The RBS and NEP survey thus comple- 
ment each other in terms of depth and angle. 

Here we combine these 2 surveys to provide new obser- 
vational constraints on the space density, as well as the X- 
ray luminosity function of non-magnetic CVs. The combined 
survey is still purely X-ray flux limited (although the flux 
limit is variable over the survey area), and yields a sample 



of 20 non-magnetic CVs. In a future paper we will consider 
magnetic CVs detected in the RBS. 

We describe the non-magnetic CV sample in Section [3] 
where we also present distance and Lx estimates for these 
20 systems. In Section[3]we describe the calculation of p and 
<!>, and in Sections U and [S] we present the results. Finally, 
the results are discussed in Section [S] and the conclusions 
listed in Section [7] 



2 THE FLUX-LIMITED CV SAMPLE 

The RBS covers half the sky down to a flux limit of 
Fx > 1-1 x 10 _12 erg cm _2 s _1 (assuming a 10 keV thermal 
bremsstrahlung spectrum), and includes 16 non-magnetic 
CVs. With the same assumed spectrum, the NEP flux limit 
varies from roughly 1.2 x 10~ 14 to 9.5 x 10~ 14 erg cm~ 2 s _1 
over the 81 sq.deg. survey area. Only 4 CVs where detected 
in the NEP, all of them non-magnetic (RX J1715, SDSS 
J1730, RX J1831, and EX Dra). The 2 surveys have no CVs 
in common, although they overlap slightly — none of the RBS 
systems are in the area covered by the NEP survey, and the 
NEP CVs are all fainter than the RBS flux limit. 

The complete sample of 20 non-magnetic CVs is pre- 
se nted in Table [p The RBS sample differs from that used 
bv lSchwopeeT al. (2002) in that we include TW Pic and ex- 
clude EI UMa. TW Pic was at first t hought to be an inter- 
medi ate polar l|Mouchet et al. 199ll; Patterson fc Mouldenl 
Il993h . but more recently [Norton et al. I (|2000T ) has argued 
convincingly that the data favour a non-magnetic nature. EI 
UMa, on the other hand , has now been shown to be mag- 
netic (|Reimer et al. I l2008l : iRamsay et al. | |200S| ; iKozhevnikovl 
l2010t ). 

We assume that none of the systems in our sam- 
ple are period bouncers (CVs that have evolved through 
the observed minimum P oro near 80 min, and that ar e 
now evolving towards longer Port.; e.g. IPaczvriskH Il98ll ). 
In the case of the long-period CVs, TT Ari, EF Tuc, RX 
J1831. 7+6511, WW Cet, V405 Peg, and EX Dra, this needs 
no justification. Eight of the short-period systems (SW 
UMa, T Leo, BZ UMa, VW Hyi, WX Hyi, SU Uma, SDSS 
J173008. 38+624754. 7, and TY PsA) have mass ratios imply- 
ing that they are not period bouncer s (as indicated by frac- 



tional superhum p period excess; see IPatterson et al. 120051 ; 
|Pattersodl2lmH 1 T For CC Scl, RX J1715.6+6856, IQ Eri, 
RBS490, TW Pic, and RBS 1411 there is no information on 
mass ratio (in fact, the orbital periods of several of these sys- 
tems are not known), and the simplest assumption is that 
they are normal, pre-period bounce CVs. 



2.1 Distance estimates 

In order to measure p and $, we need distance estimates 
for all CVs in the sample. Although good distance measure- 
ments are available for some of these systems, in many cases 



1 It has been suggested that VW Hyi has a substellar donor 
star (on the basis of its near-I R spectrum; see [ Mcnnickc nt et all 
l2004h ; however, the mas s ratio dPattersorJll998| ) and white dwarf 
mass llSmith et alj|2006l) do not support this, and Hamilto n et al.l 
ll201ll) show that the donor spectral type is no later than M9V. 
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the distance estimates that we can derive are quite impre- 
cise. 

T Leo, BZ UMa, RBS490, SW UMa, and V 405 Peg have 
high- q uality parallax distance measurements jTh orstcnscn 



20031; IThorstensen et all 120061 ; iThorstensen et all 120081; 
Thorstensen et alJl2009F 7 



Other reliable estimates of CV distances are found by 
photometric parallax, in cases where the white dwarf or the 
donor star is detected in a way that makes it possible to 
disentangle its flux contribution from that of other light 
sources. We use such estimates for TT Ari and EX Drfl 
TT Ari has a distance measuremen t based on the donor 
star spectrum (Ga nsicke et al.lll999l ; this is also consistent 
with the distance derived from the white dwarf spectrum). 
Based on the photometric parallax of the secondary, the dis- 
tance to EX Dra is 240^^ pc (seelShafter fe Holland! 120031; 
iBaptista et al.ll2000l ; iPretorius et al.ll2007bl ; we have revised 
this to be slightly less conservative than the estimate we 
used previously). 

There are a few well-known, less direct (and less reli- 
able) methods of estimating distances to CVs. We u se meth- 
ods b ased on dwarf nova (DN) out bur st maximum ([Warner] 
Il987l ; see also lHarrison et al.|[20o3 and lPattersoifeOllI ). the 
near-IR apparent brightness for systems in which the donor 
star is not directl y detected (th e m ethod of Bailey 198 j . but 
as prescribed by iKniggd 120061 and lKnigge et al.||201ll) . and 
the strength of H/3 emission lines (|Pattersonl 1 19841 ; see also 
IPattersonl [20111 ). 

The relation between the absolute magnitude at DN 
outburst maxim um and P or b was most recently studied by 
IPattersonl l|201lf) . who confirms that the scatter is relatively 
small, and that there are no large outliers. We therefore use 
this relation as far as possible (for CC Scl, VW Hyi, WX Hyi, 
SU UMa, TY PsA, WW Cet, SDSS J1730, and EF Tuc). It 
should be noted that SU UMa has a parallax distance esti- 
mate that agrees very well with the distance based on out- 
burst maximum, but which is less precise; we choose in this 
case to use t he estimate based o n ou tburst, rather th an the 
parallax ( see IThorstensen! 120031 and IPattersonl |20 1 if ) . Two 
systems, SDSS J1730 and EF Tuc, do not have orbital in- 
clination measurements (although it is known that they are 
not eclipsing), and have less well determined maximum ap- 
parent magnitudes, leading to more imprecise distance esti- 
mates from this method. 

For RX J1 831 we use the pres cription of lKniggel (|2006t ) 
(as updated bv lKnigge et al.ll201ll ). This is based on a semi- 
empirical donor sequence for CVs, and the offsets between 
this sequence and the absolute JHK magnitudes of a sam- 
ple of CVs with parallax distances. Apparent near-IR mag- 
nitudes for RX J1831 we re obtained from the Two Micron 
All Sky Survey (2MASS; Iskrutskie et al.ll2006l ). 

Finally, although there clearly exists an empirical rela- 
tion between EW(H(3) and the absolute m agnitude of the 
disc, this relation contains large scatter (see IPattersonl l201ll 
for an updated plot). We therefore use it as a last resort, in 



2 ISion et"af I 1119951) estimate a distance of 75 pc for VW Hyi, 
based on FUV spectroscopy. However, since it is not clear what 
contributi on sources other tha n the white dwarf make to the FUV 
light (see iGodon et al.ll20oi) . we prefer to disregard this value 
(although it is completely consistent with the estimate we will 
use). 



those four cases where the data required by the other two 
methods are not available (TW Pic, IQ Eri, RX J1715, and 
RBS1411). For RX J1715 and RBS1411, we have no data to 
allow us to check whether the absolute magnitudes we find 
are reasonable (other than that RX J1715 is known to be a 
short-period CV, and that RBS1411 has an optical spectrum 
resembling that of a short-period CV). For the remaining 2 
systems, we find absolute magnitudes that are in reasonable 
agreement with what we would find from outburst (in the 
case of IQ Eri), and the method of Knigge (in the case of 
TW Pic), if we made reasonable assumptions about orbital 
periocjf|. 

Interstellar extinction is expected to be low for our sys- 
tems, since they are at high Galactic latitude, and since most 
of them are quite nearby. For those systems with distances 
below 200 pc, we neglect interstellar extinctio n. For the 
more distant objects, we use Ay estimates from IPattersonl 
i|201ll ), where available, and f or a few more, we find esti- 
mates in lBruch fc En gel (1994). In the cases where no more 
di rect estimate of ex t inctio n is available, we use the model 
of|A mores fc Lepinel ([20051 ) . with a few iterations, so that 
the value we finally adopt in the distance calculation is that 
given by the model at the estimated distance to the ob- 
ject. To convert from visual extinction to extinction in the 
2MASS bands, w e use Aj = 0.282^4^, A H = Q.17bA v , and 
A Ks = 0. 112,4 v l|Cambresv et al.i r2002). We conservatively 
assume errors of 50% in the extinction values. 

We then find the probability distribution for the dis- 
tance to each source, assuming Gaussian errors in all 
the input parameters (apparent magnitudes, inclination, 
EW(Hf3), extinction). The distances estimates, listed in col- 
umn 3 of Table [1] are in all cases the median, together with 
the 1-ct confidence interval corresponding to the 16th and 
84th percentile points. 



2.1.1 Possible bias in the distance estimates 

Distance estimates may suffer from t wo well-known biases, 
namely Malmquist (|Malmquistiri924l ) and Lutz-Kelker bias 
l|Lutz fc Kelkerl 1 19731 ). These biases, the relation between 
them, and how to c orrect for them have been the subjec t 
of many papers (e.g. iGonzalez fc Faber|[l997l ; ISmithll2003T ). 
Here we will examine whether our distance estimates could 
be biased. 



3 IQ Eri has been observed in dwarf nova outburst; the quiescent 
spectrum, large outburst amplitude, and the apparent long out- 
burst interval all suggest that this is a short-period system, so it 
is possible to find a very rough indication of distance based on 
the outburst magnitude, which is in satisfactory agreement with 
the distance based on EW(H/3). TW Pic has spectroscopic and 
photometric periods near 2 and 6 hours (the orbit al period was 
initially assumed to correspond to the 6-h signal; iNorton et al.1 
(2000) discuss the alternative interpretation of the more stable 
2-h period being t he orbital peri od). Assuming the period near 2 
hours is P or b, the iKniggel ||2006l ) sequence would give a distance 
in agreement with what we find from EW (H /3) . However, if the 
orbital period is as long as 6 h, the observed K = 14.1 would 
imply an improbably bright absolute V magnitude (My ~ 3.2; 
although note that at such long P or b, the system would proba- 
bly have an evolved donor, for which the donor sequence is not 
appropriate). 
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Table 1. The 20 non-magnetic CVs detected in the RBS and NEP survey, together with their orbital periods, distances, and 
X-ray fluxes and luminosities. The meaning of 1 /Vj and pj / po is explained in Section l3,2.1l References are for published distances, 
or the values of P or t, EW(Hj3), and binary inclination used in estimating distances. 
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Lutz-Kelker bias affects parallax me a surem ents. This 
was carefully considered by IThorstensen! {2003), and the 
same procedure was used by Thorstensen (2006, 2008, 2009). 
We are therefore confident that the parallax distance esti- 
mates used for 5 of the CVs in our sample are unbiased. 
This leaves the possibility of Malmquist bias in the distance 
estimates of the remaining systems. If some type of objects 
have average absolute magnitude (M), with intrinsic scatter 
(Jm, and the apparent magnitude of such an object is used to 
estimate its distance, the distance is biased for a magnitude- 
limited sample, because the absolute magnitude distribution 
of the magnitude-limited sample has a mean brighter than 
<M). 

We will assume that the calibrations we use for the dis- 
tance estimates based on DN outburst maximum, EW(Hf3), 
and near-IR apparent brightness are not biased. This is rea- 
sonable, first because the uncertainties are large compared 
to any expected bias, but also because the samples used 
to derive these relat ions are not s eriously affected by flux 
limits. For example, iKni ggc ( 2006)) used CVs with reliable 
parallax distances to determine the offset between his donor 
sequence and the absolute magnitudes of the CVs, and only 
one out of the 23 systems was excluded by the flux limit 
of 2MASS. The same is true for the other two relations — 



they are derived from samples of CVs with well-measured 
distances, and these samples are not strongly influenced by 
a flux limit. 

An important point is that our distance estimates are 
based on optical apparent magnitudes (or, in one case, near- 
IR apparent magnitude), while the surveys are X-ray flux- 
limited. It is easy to verify that, if optical- and X-ray lumi- 
nosity are uncorrelated, there is no bias the distance esti- 
mates. The bias only appears if L op t (or Lie) is correlated 
with Lx- A simple expectation is that L opt is proportional 
to Lx, because both should be proportional to M; however, 
this is not true. The rel ation between Lx and L op t is time- 



dependent for DNe (e.g. Jones fc W atson 1992), and flattens 
off at bright L opt (e.g. IPatterson fc Ravmondll 1985a! ) . These 
complications aside, the bias will be strongest for Lx oc L opt ; 
we therefore use this assumption in determining how seri- 
ously our distance estimates might be affected. 

The strength of the Malmquist bias clearly depends on 
the intrinsic scatter in absolute magnitude, <tm, but also 
on several factors that combine to determine how close the 
flux-limited sample is to a volume limited sample; these are 
6, Galactic scale height h, absolute magnitude M, distance, 
and the survey flux limit. We consider each of the 15 CVs 
that may have distance estimates suffering from Malmquist 
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Table 2. The 5 CVs with distances that might be biased. We give 
the estimated distance, and the factor by which the real distance 
might exceed this if optical- and X-ray luminosities are correlated. 



multi-temperature (or cooling flow) models a re needed to 
provide satisfactory fits to observations (e.g. Mukai et alj 



System 


draw 1 PC 


^corrected/ ' &raw 


CC Scl 
EF Tuc 
TW Pic 
RBS1411 
IQ Eri 


359 133 
346+iS 

U6±£« 


1.26 
1.20 
1.55 
1.51 
1.72 



bias in turn, using the appropriate b, {M), and am, together 
with the X-ray flux limit of the survey it was detected in, 
and an exponential vertical density profile for the GalaxjQ. 
Using a Monte Carlo simulation for each system, we itera- 
tively find the distance at which a population with appropri- 
ate Gaussian M distribution, subject to an X-ray flux limit, 
would result in an estimated distance distribution with me- 
dian equal to our distance estimate. 

As expected, the bias is not present for sufficiently small 
distance or gm, or sufficiently deep X-ray flux limit. With 
the assumed correlation between L opt and Lx , we would es- 
timate biased distances for only 5 of the CVs in our sample 
(listed in Table [5}. Although a few of the distance estimates 
may be seriously biased, these are systems that make only 
very small contributions to the space density and luminosity 
function (see the 6th and 7th columns of Table [T] and the ex- 
planation in Section[3]2} ■ Correcting these 5 distances for the 
possible bias does not significantly change the results that 
will be presented in Section 2] and [S] Therefore, although we 
do not know the relation between Lx and L opt , and thus 
whether in principle we should correct the distances, this 
possible bias can be safely neglected. 

However, even if unbiased distance estimates are used, 
the resulti ng luminosity fun ction still contains Malmquist- 
type bias (|Stobie et alj ' ll989h . 
tionIO 



We will return to this in Sec- 



2.2 X-ray luminosities 

Most non-magnetic CVs have X-ray emission that can be de- 
scribed as thermal bremsstrahlung from the boundary layer 
(the inner part of the disc where the flow is no longer ke- 
plerian, but is slowing down to match th e rotation of the 
white dwa rf; Patter son fe Ravmondlll985al . but see also e.g. 
iPerna et al]|2003l ). The X-ray spectrum should then be the 
sum of emission from material with temperatures ranging 
from the temperature of the shock at the outer edge of 
the boundary lay er to the tempera ture of the white dwarf 
photosphere (e.g. iMukai et al1l2003f ). X-ray observations of 
many systems can be fit by a single temperature ther- 
mal bremsstrah lung spectrum with kT between roughly 5 
and 20 keV fe.g.lPatt crson fc Raymond! 1985al ; I Vrtilek et al 
1994 lMukai_etaljll997l; ISzkodv et~afl IgOOOl: [Baskill et al 



20051 ; Pandcl ct al. 2005; Mu kai et alj|2009h . In other cases 



4 Our simple Galaxy model is discussed further in Section 13.11 
Here we use a scale height of 260 pc; this is perhaps too large for 
younger, long-period CVs, but it is a conservative assumption, 
since the strength of the bias increases with scale height. 



Baskill et alTEoO 
Hoard et al] 1201 



| 2003l; iPandel Cordova fc Howelll |2003| : 
Pandel et alj 120051 iHilton et al.l [ 20071 . 
iBvckling et al.l l201Ch . High-M systems (nova-like CVs in 
a high state and DNe in outburst) are expected to have 
optically thick boundary layers, and therefore much softer 
spectra; however, part of the boundary layer remains op- 
tically thin so that the soft compo nent does not dominate 
the spectrum above ~ 0.5 keV (e.g. IPatterson fc Raymond! 



1985al; IPatterson fc Raymond! f985bl ; |Jones fc Watsonlll99a 
Wheatlev et alj|2003t ) 

We assume a kT = lOkeV thermal bremsstrahlung 
spectrum for all CVs in our sample, and quote X-ray fluxes 
and luminosities in the 0.5-2.0 keV band. Although there 
is no good physical justification for this simple approach, 
it is acceptable for our purposes, because the energy band 
we are using is narrow (decreasing the sensitivity of Lx to 
the assumed spectrum), and because our distances are for 
the most part quite imprecise (implying that error arising 
from the assumed spectral shape is unlikely to contribute 
significantly to the to tal error in Lx', see b elow) . We use the 
dust to gas ratio of IPredehl fc Schmittl (|l995h to convert 
the Av estimates to Nh for each system. The Nh estimates 
are low (the hig hest being 5.4 x 10 20 cm -2 , for SDSS J1730; 
IPattersonfeoill i. We list unabsorbed F x and L x in Tabled] 
In order to provided further justification for our as- 
sumptions regarding the X-ray spectrum, we can consider 
the error in Lx using standard error propagation: 



~L X 



<7F X 



F x 



Where ad is the error in distance (in the best cases about 
15%, but larger for most systems), of x obs is the obser- 
vational error (typically ~10%), and op x c is the error 
caused by uncertainty in Nh, and by an incorrect assump- 
tion of X-ray spectrum. In order for of x spcc to dominate 
ol x , we must have op x s cc /Fx > 0.3. This seems unlikely, 
because in the narrow energy band (0.5-2.0 keV), for mod- 
erate Nh and allowing for a 50% error in Nh, the difference 
in Fx between single-temperature bremsstrahlung spectra 
with kT = 2keV and 20keV is only ~ 20%. 

High-M CVs are naturally expected to have brighter 
Lx', therefore, Lx should increase with P or b- This is cer- 
tainly not apparent in the l uminosities we find (see t he P or b 
and Lx values in Table [TJ. IPatterson fc Raymond! (|l985ah 
show that, for low-M systems, Lx indeed increases with M, 
but that at roughly 10 16 g/s the relation between Lx and M 
flattens off. Binary in clination may also b e expected to add 
noise to this relation. iBaskill et al.l (|2005 | ) find only a weak 
correl ation between Lx and P or b, while Ivan Teeseling et all 
( 1996) find none. Given the large scatter or weakness of any 
correlation between Lx and P or b, Lx is not a good indica- 
tor of period (or, by implication, age); in the calculations 
described in Section 15.21 Lx is therefore the only physical 
parameter of the simulated CVs. 
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3 MEASURING THE SPACE DENSITY AND 
LUMINOSITY FUNCTION 

Here we describe the calculation of the space density and X- 
ray luminosity function from the sample of 20 non-magnetic 
CVs detected in the combined RBS and NEP survey. The 
results of this calculation are presented in the next section. 

3.1 Approximations and assumptions 

The most important assumption of the method we use is 
that the observed CV sample is representative of the under- 
lying population. Whether this is a reasonable assumption 
is addressed further in Section [5] and [6] 

The number of CVs per unit volume is of course a func- 
tion of position in the Galaxy. We ignore the (weak) radial 
dependence of p, and assume that the vertical density profile 
is exponential 

\z\/h 



p(z) = p e" 



(1) 



with z the distance from the Galactic plane (i.e. z — dsinfe). 
The local space density is denned as po = p(0), the mid- 
plane value of p. 

The X-ray luminosity function, $ (logLx) is defined so 

that 



Po,l x = $ (logLx) dlogLx, 



(2) 



where po,L x is the local space density of CVs in the luminos- 
ity bin of width dlogLx centred on logLx- In other words, 
po is the integral of $ over logLx • 

In calculating po, we will compare two assumptions re- 
jar ding Galactic scale height. First, as in iPretorius et al.l 
2007bf i. we take h = 120 pc for long-period CVs, and 260 pc 
for short-period systems, as a (quite crude) approximation 
of the fact that these are young and old populations, re- 
spective!^. The second po calculation uses a single scale 
height of 260 pc for all systems, and gives results that are 
not significantly different (see Section 271]). This implies that 
the simpler approach of using only a single scale height is 
justified, and that is what we will assume in calculating $ 
(the reason we require the simpler assumption for $ will be 
explained in Section \5. 21) . 

Interstellar extinction is computed by integrating the 
density of the interstellar medium along each line of sight to 
obtain the neutral hydrogen column density, Nh- We again 
assume only vertical dependence in the density of gas; i.e., 

Pism(z) = pi SMfi e-^ /hlSM . (3) 

We take h ISM = 140 pc (e.g. iRobin et all 120031 ; 

iDrimmel fc Spergell [2001 1. and we find Pis mo by assum- 
ing A v = 2mag/kpc for b = deg dAllenl ll976Tpl . and 
N H = 1.79 x 10 21 cm- 2 A v l|Predehl fc Schmitt) n"995V 



5 This is expected theoretically (e.g. iKolb fc Stehl j Il996l) . but 
observations have not yet clearly shown a differe nce in scale height 
between long- and short-period CVs (see e.g. ISzkodv fc H owell 
1 19921 ; Ivan Paradiis et al.lll996l ; iNorth et al1l2002l ; lAk et al.ll20ld) . 

6 Some more recen t work indica t es low er midplane extinc- 
tion values; e.g., IVergelv et_ah I lll99Sh fin d 1.2m ag/kpc, 
while the models of [ A mores fc Lepinel j2005f) and 
IDrimmel. Cabrera-Lavers. fc Lope^Corredoiral i2003h give 
lower average values still. However, the results that will be 
presented in the next two sections are not sensitive to this 



As discussed in Section 12.21 we will assume a kT = 
10 keV thermal bremsstrahlung spectrum for all CVs in our 
sample. We will explain below that the p calculation uses the 
maximum distance at which a given CV could have been de- 
tected. Since this maximum distance depends on the ratio 
of Fx to the flux limit, which may as well be expressed as a 
ratio of count rate to limiting countrate, the assumed X-ray 
spectral shape has very little influence on the p calculation 
(only the interstellar absorption depends on the X-ray spec- 
trum). 



3.2 The calculation 

Since the combined survey is complete (up to the variable 
flux limit), we simply need to count the systems detected 
inside the observed volume to calculate p (and do the same in 
a set of luminosity bins to find <!>). However, the dependence 
of p on z, as well as the flux- rather than volume-limited 
nature of the sample need to be accounted for. Both these 
complications are taken care of by the well-known 1/Vmax 



method (see e.g 



ll989l : lTinnev et alJll993fi 



Schmidt! [l96g|; iFelterJ Il976l; IStobie et all 



Our aim is to find not only a best estimate of po and 
but also to estimate the errors on both, accounting for all 
sources of statistical error (including distance errors, small 
number statistics, and observational error on Fx or count 
rate). We describe in the following sections first how we 
implement the 1/V ma x method, and then the Monte Carlo 
simulation to find the uncertainties on our measurements, 
as well as the numerical tests we use to verify that the error 
estimate for po is correct. The uncertainty in $ proves to 
be harder to estimate, and we will return to this issue in 
Sections rOl and [531 



3.2.1 Basics of the 1/V ma x method 

The 1/V m ax method essentially allows the 'volume limit' 
of the survey to vary according to the luminosities of the 
systems in the sample. For each observed system, we find 
the 'generalized' (i.e., taking into account the exponential 
vertical density profile) maximum volume, 



y. = n fe3 

3 jsin&| 3 



[2- (x] + 2x ] +2) fT*i] 



(4) 



(e.g. iTinnev et aDll993h . The index j represents the CVs in 
our sample; f2 is the solid angle covered by the survey and 
Xj = dj\sinb\/h, with dj the maximum distance at which 
CV j could have been detected (given its luminosity and 
the survey flux limit) . Equation [3] assumes the spatial de- 
pendence of p as given by equation [1] Vj is then the volume 
probed by the survey for sources with Lx the same as the 
observed system j. Because b (and in the case of the NEP 
survey, also the flux limit) is variable over fl, we compute 
each Vj as a sum over smaller solid angles, 8£l. For the RBS, 
each 5£l is a slice subtended by 2° in b (excluding the small 
area that is also covered by the NEP survey) , while the NEP 



assumption, and do not change significantly if we take the 
midplane extinction to be as low as Ay = 0.75mag/kpc. This 
is probably because the two surveys are restricted to quite high 
Galactic latitudes. 
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area is divided into a 36 x 36 pixel grid. For both surveys, 
the 80, are sufficiently small that the error introduced by b 
varying over Sil is negligible. The space density po is then 
the sum of the space densities represented by each CV, i.e. 

The luminosity function $ in a given luminosity bin is cal- 
culated similarly, by restricting the sum to systems in that 
luminosity bin. 

The 1/Vj and Pj/po values for the systems in our sam- 
ple, obtained with the maximum distance found by using the 
best d and Fx estimates, are given in Table Q] These values 
thus neglect for the moment all uncertainties. We list them 
only as an indication of the contribution that each system 
makes to the total po. 



3.2.2 The Monte Carlo simulation 

In order to find the error on po, we compute its probabil- 
ity distribution function using a Monte Carlo simulation 
that finds po (as described above) for a large number of 
mock samples with properties that fairly sample the param- 
eter space allowed by the data. This also produces a best- 
estimate The appropriate mock s amples are generated as 
described in lPretorius et al.l ( |2007bh . 

Briefly, we treat each observed CV as a representative 
of a population of similar systems, and allow each popula- 
tion to contribute to each mock sample. For a given mock 
sample, we generate 20 mock CVs by drawing Fx and d 
for a mock CV from the probability distribution functions 
(Gaussian errors in Fx and d distributions as computed in 
Section [2. II) of a real observed system; thus each mock CV is 
the counterpart of one of the observed CVs. This takes care 
of the uncertainties in Fx and d (or, combining those two 
parameters, Lx). The Poisson uncertainty associated with 
the small sample size is accounted for by a weighting factor, 
p, drawn from the probability distribution of the number of 
sources belonging to the population (corresponding to a par- 
ticular observed system) that one expects to detect in the 
combined survey. Using Bayes' theorem, this distribution is 
related to the Poisson distribution P(N t, s \p) by 

P(p\N obs ) cx P(N obs \p)P(p), 

where the actual number observed in the RBS and NEP sur- 
vey is Nobs = 1, since we take each observed CV as represent- 
ing a p opulation of similar systems. As in IPretorius et al.l 
(2007b), we use the uninformative prior P(p) = 1/p. Then 
for each mock sample, po = ^2jPj/Vj, and $(logLx) is 
found by summing only over the j corresponding to mock 
CVs with luminosities placing them in that logix bin. We 
generate a large number of mock samples, and calculate pa 
and $ for each of them. We then add the $ values from all 
the mock samples in a given bin, and divide by the number 
of mock samples. The distribution of po is normalized to give 
a probability distribution function. 

In our po calculation, the largest part of the error budget 
comes from uncertainty in distance (or, equivalently, Lx). 
The calculation naturally yields the error on po, but the 
situation for $ is more complicated (see Section \3. 2. 30 . 



3.2.3 Numerical tests 

We carried out several tests to confirm that the Bayesian 
procedure outlined above produces reliable error estimates 
for po- The input for each test is a Galaxy model (the same 
as assumed by equation [T] and survey area and flux limit, 
and a CV luminosity function. We experimented with sev- 
eral different luminosity functions, as well as with including 
simulated errors in luminosity for the "observed" CVs. 

A given test first simulates an "observed" sample, and 
then calculates the po distribution based on that sample 
in the same way as described above. This is repeated to 
generate many samples with their po distributions (note that 
for a given input, the samples do not all have the same 
size). For certain input luminosity functions, we found that 
for ~68% of these simulated samples, the true input po is 
contained a 1-a interval of the calculated po distribution. 
This holds when the combination of input flux limit and 
luminosity function produces an average of 20 sources per 
simulated sample, but also when the simulated samples are 
much larger, and the po distributions become narrower. 

The luminosity functions for which the tests fail (in 
other words, for which we cannot recover the input po) are 
those in which a significant contribution to p is made up by 
systems so faint that none of them are included in the "de- 
tected" samples. This is then what it means for the observed 
sample to be representative of the intrinsic population — at 
least one CV belonging to any intrinsically large population 
is detected. 

We perform similar tests to see if we can recover the lu- 
minosity function used as input. The result is that we can, 
but only when the errors in luminosity we assign to "ob- 
served" CVs are small compared to the size of the luminos- 
ity bins. In this case, we can always reliably place a given 
CV in a single Lx bin, and the error on each <E> bin can be 
found in the same way as the error on po. However, clearly 
also here, we cannot recover any faint end of the luminosity 
function beyond the faintest "detected" Lx- 

When we simulate a survey where the errors on Lx of 
the "observed" CVs are large, we recover a luminosity func- 
tion that is broader and flatter than the input. Furthermore, 
we can no longer find the error on a given $ bin as before, 
from the distribution of values in that bin. The reason for 
this is the correlation between pj and Lx for that system. 
A fainter system can be detected only to a smaller distance, 
and hence yields a larger 1/V); this is illustrated in Fig. [T] 
To understand why the calculation does not give the error 
on $, one can consider constructing a luminosity function 
with bins sufficiently small that, for each mock sample, ei- 
ther 1 or systems fall in a given Lx bin. For some bin of 
this luminosity function, the value of <&, if it is not 0, will 
be determined by the correlation shown in Fig. [T] with rel- 
atively small scatter, determined by the weighting factor p. 
The most important source of error will be gone, since Lx 
is known, in a given bin. 

For our real sample, the uncertainties in Lx are large 
enough that a given mock system is not always in the same 
Lx bin (even for very coarse bins), and in part because 
we compute $ over a large range in logLx, a typical mock 
sample leads to many empty bins. This implies, in a given 
bin, many points (corresponding to many mock samples) at 
$ = 0. The distribution of the non-zero values is determined 
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Figure 1. The correlation between pj and Lx for the RBS (fine 
curve), the NEP survey (dotted curve), and the combined survey 
(bold curve). This shows the pj one would find for a given system, 
as a function of its Lx ■ The correlation becomes flat at the bright 
Lx end, because a sufficiently bright system can be detected out 
to the edge of the Galaxy. 

by the correlation, with perhaps a factor of 2 or so spread (1 
or 2 systems in that bin) and with some additional scatter 
from the distribution of the weights. Therefore, the distri- 
bution of values in a given bin cannot tell us the error on $ 
in that bin, and more work is needed (see Section [5. 2[1 . 



4 ESTIMATES OF THE SPACE DENSITY 
AND X-RAY LUMINOSITY FUNCTION 

Here we present the results of the calculations described in 
Section[3]above. We will elaborate on the interpretation and 
limitations of these results in Section [S] 



4.1 The probability distribution function of po 

The probability distribution function of po resulting from 
the calculation described above is shown in Fig. As al- 
ready noted, we find that the distances errors dominate 
the total uncertainty in po, although the small sample 
size also contributes significantly. The mode, median, and 
mean of the distribution are 2.3 x 10~ 6 , 4.4 x 10" 6 , and 
7.4 x 10 -6 pc~ 3 , and are marked by solid lines. The dashed 
lines at 2.2 x 10~ 6 and 1.0 x 10~ 5 pc -3 show a 1-a confi- 
dence interval (the 16th and 84th percentile points of the 
distribution). 

Considering the large errors, our p estimate of 4.4^2 x 
10 -6 pc -3 resulting from the combined RBS and NEP sur- 
vey is in reas o nable agreement with the value found by 
ISchwope et al.l (|2002T ) when omitting the two systems for 
which they estimated very small distances. We will compare 
the new estimate with the higher value we found from the 
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Figure 2. The po distribution resulting from our simulation. For 
this distribution, we used scale heights of 260 pc and 120 pc for 
short- and long-period systems, respectively. Solid lines mark the 
mode, median, and mean at 2.3 X 10 -6 , 4.4 X 10 -6 , and 7.4 X 
10 _6 pc — 3 . Dashed lines show a 1-a interval from 2.2 X 10 — 6 to 
1.0 X 10~ 5 pc -3 . 

NEP survey alone in some detail in Section 15.31 Consider- 
ing short- and long-period CVs separately (even assuming 
that TW Pic, RBS490, RBS1411, and IQ Eri are all short- 
period systems), we find that the space density estimate is 
dominated by long-period systems; po of long-period CVs 
is roughly 1.2 times that of short-period systems (this is 
discussed further in Section [6]). 

If we use a single Galactic scale-height of 260 pc for 
all systems, we find po = 3.8 jl^'l x 10~ 6 pc -3 ; this is not 
significantly different from the result when using different 
scale heights for long- and short-period CVs. Similarly, if 
we correct for the possible distance biases listed in Table f2J 
the resulting po is insignificantly different, namely TOIj'q x 
10 _6 pc" 3 . 

4.2 The observed luminosity function 

The same calculation that gives po also yields as explained 
in Section 13.2.11 The X-ray luminosity function obtained 
from this calculation (assuming h = 260 pc for all CVs) is 
shown in Fig. [3] We also plot there the log(Lx/erg s _1 ) es- 
timates of the 20 detected CVs, in order to show the size 
of the errors in these values. The uncertainty in Lx has the 
effect of smoothing features in the luminosity function. This 
explains, for example, why the inferred luminosity function 
has a gradual turn-over at the faint Lx end, rather than 
a sharp cut-off at the faintest observed Lx value. If we in- 
tegrate this histogram, we find 6.1 x 10 _6 pc -3 , a number 
between the median and mean of the po distribution. 

As explained in Section f3. 2. 31 this calculation does not 
provide the uncertainty on hence we do not show er- 
ror bars on the histogram in Fig. [3] Furthermore, although 
the result does not change significantly if we correct for the 
possible distance biases discussed in Section ^. 1.11 the lumi- 
nosity function probably still suffers from Malmquist-type 
biases, and while the turnover at the faint end might rep- 
resent a real cutoff, it might also just reflect the limited 
depth of the observations. We will return to these problems 
in Section \5. 2 1 



Space density and luminosity function of CVs 9 




29 30 31 32 

log(L x /erg s" 1 ) 



Figure 3. The histogram shows the observed X-ray luminosity 
function. As explained in Section 14.21 the calculation that yields 
this histogram does not produce an estimate of the errors on it. 
The points with error bars show the estimates of Lx for the 20 
non-magnetic CVs that the calculation is based on (these points 
are arbitrarily offset in the vertical direction, for display pur- 
poses). 

5 THE EFFECT OF THE FLUX LIMIT 

Whether a given survey produces a CV sample that is rep- 
resentative of the underlying population is determined by 
its depth and area. No evidence of even a large population 
of CVs is expected to show up in a a flux-limited survey, 
if such a population is sufficiently fainfl We can, however, 
still place a limit on the size of a hypothetical faint popu- 
lation of CVs that may have gone undetected in these two 
surveys; this is done in Section \5. li b elow . In detail, the effect 
of the flux limit on the p estimate that a given survey pro- 
duces depends on how steeply the luminosity function rises 
towards the faint end and on where it cuts off or drops. In 
Section [5.21 we return to the matter of errors on <E>, and at- 
tempt to constrain the shape of the true luminosity function 
that, after being subjected to the appropriate observational 
biases and errors, would give rise to the (smoothed) lumi- 
nosity function we actually observe. In Section [5.31 we also 
compare the two individual surveys that were combined to 
construct our sample, to verify that they are consistent with 
each other. 

5.1 Upper limits on the space density of an 
undetected population 

In order to calculate how large a faint CV population could 
have gone undetected, we use another Monte Carlo simula- 
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Figure 4. The upper limit on the mid-plane space density as a 
function of X-ray luminosity for an undetected population of CVs. 
The data from the simulation are shown as a bold histogram, and 
a fit is over-plotted as a finer line. The 2 upper, fine histograms 
show the corresponding results for the RBS (middle) and NEP 
(top) surveys alone (the top histogram is what was presented as 
figure 7 of Pretorius et al. 2007b, except that we here take the 
Galactic scale height as 260 pc, in order to be consistent with 
what we do for the RBS survey and the combined survey). 

tion. We assume that there is a population of CVs, all with 
the same Lx, and find the upper limit on its space density, 
based on detecting no member of this population. While the 
assumption of a singie-Lx population is unphysical, it gives 
useful results that are easily expressed in terms of Lx- 

We again model the Galaxy as having an exponential 
vertical density profile and no radial dependence in the num- 
ber density of stars; we also include extinction, in the same 
way as described in Section 13.11 We assume a single scale 
height of 260 pc for the hidden population (as we did in 
Section [4.21 when calculating $). We then find the value of 
po for which the predicted number of detected systems is 3 
(detecting systems such systems is then a 2-cr result). 

Fig. [4] shows the maximum allowed po as a function of 
Lx for CVs that make up the hypothetical hidden popula- 
tion, from the RBS and NEP survey separately (the middle, 
and top-most fine histograms, respectively), as well as for the 
combination of the two surveys. The limit from the simula- 
tion for the combined survey is plotted as a bold histogram, 
and the fine curve is a fit to the data, given by 

Pmax = 4.91 x lO-^Lx/lO^ergs- 1 )- 1 - 39 ?^ 3 . 

Thus for = 3.6 xlO 28 ergs" 1 , we have p < 2xl0" 4 pc" 3 . 
This limit is stronger than we previously found for the 
NEP survey al one (po of 2 x 10~ 4 p c~ 3 implying Lx < 
2 x 10 29 erg s~ 1 ; IPretorius et al.ll2007rl ; this is a result of the 
much larger volume probed by the RBS (see Section |5.3[) . 
Note that here po refers only to a possible undetected popu- 
lation, and does not include the small contribution from the 
observed systems. 



7 A volume limited survey, on the other hand, might completely 
miss CVs that are intrinsically bright and rare, if the volume is 
small; however, rare systems are by definition not very important 
in finding the size of the underlying population. 



5.2 Forward modelling to account for bias in $ 

As noted earlier, the luminosity function presented in Sec- 
tion 14.21 is expected to suffer from Malmquist-type biases, 
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and we were not able to estimate the uncertainty on it. 
Biases in luminosity functions comp uted by the 1/Vma x 
method are discus sed in detail by e.g . IStobie et al.l (11 989). 
iGeiio et all |2006l ), and iTorres et all |2007l ). IStobie et al l 
( 19891 ) show that even when unbiased distance estimates are 
used, $ is still biased. In order to correct for these effects, 
we have to assume a functional form for the true luminos- 
ity function. We will model it as a power law, and then 
determine the power law index that best reproduces the ob- 
servations we have. 



5.2.1 Monte Carlo simulation of the survey 



In the same way as described in Section f3. 2. 3 1 we populate a 
model galaxy with CVs from an input luminosity function; 
here the form is 

<t = const x L x 

for the range 28.2 < log(L x /erg s" 1 ) < 32.0, and $ = 
elsewhere (in practice, of course the input $ is discrete). 
Using the flux limit and area of the combined RBS and NEP 
survey, we then find the sample that would be detected and 
calculate an output $ exactly as in Section ^. 21 Errors in Lx 
similar to those of the real CV sample are assigned to the 
"detected" CVs (these errors are Gaussian in log(Lx), which 
is not quite the case for all the real systems; see Table [T]). 
For a given input "I>, we repeat this many times (note that 
the "detected" CV samples do not all have the same size). 
Then we compare the distribution of the output in every Lx 
bin with the observed $ presented in Section \4. 21 

We vary the power-law index a, to find the best input 
$. The normalization constant is fixed (for a given a) so that 
the integral of $ over the range 29.8 < log(Lx/erg s _1 ) < 
31.8 is the same for the assumed input as for the observed 
$, plotted in Fig. [3] The position of the faint cutoff in the 
input $ is not important — we would obtain the same result 
if it extended to fainter Lx- However, allowing the input "!> 
to be non-zero to log(Lx/erg s _1 ) > 32 results in output 
values much higher than the data in the brightest few Lx 
bins (the brightest position of the cutoff consistent with the 
data depends on the the power law index, but we performed 
only a single-parameter fit, with both bright and faint cut- 
off Lx fixed, because the calculations are computationally 
expensive). We use a single scale height in these simulations 
(260 pc), because Lx is the only property that a simulated 
CV has (we cannot assign an age, because, as mentioned 
before in Section [2.21 Lx is not a clean indicator of P or b). 

Since the "detected" CV samples in these simulations 
are subject to the same selection criteria as the real sample, 
they are affected by Malmquist-type biases in the same way 
as the data. Therefore, to the extent that our simple Galaxy 
model is suitable, the best-estimate power-law index, found 
by comparing these simulations to the observations, is un- 
biased. 



5.2.2 Best- estimate power law $ 

The assumed <3> resulting in the best-fit model output is 
shown in Fig. O together with the output from the sim- 
ulation, as well the observed $ from Section 14.21 We find 




30 31 
log^/erg s" 1 ) 

Figure 5. The assumed intrinsic $ (crosses), and the output 
resulting from the calculation described in Section 15.2.11 (points 
with error bars), plotted over the observed <I> (histogram). The 
data shown here are the input that gives the best match between 
the simulated output and the observed The points and error 
bars are the median and a 1-cr interval of the the distribution of 
values in each bin, resulting from many "detected" samples. 



a = —0.80 ± 0.05, where the error on a is based on \~ m ~ 
creasing by 1. For this best-fit, the reduced \ 2 is 

Note that the output is lower than input <E> at the faint 
end; this is the expected effect of the flux limit. Also, at the 
bright end, the sharp cutoff in the input is not recovered, 
because of error in Lx- We expect that the errors on the 
output from these simulations give a reliable indication of 
the uncertainty on the observed 



5.2.3 The distribution in z 

This calculation also allows us to check the distribution in 
height above the Galactic plane of the CV sample. Since 
it is a high Galactic latitude, flux-limited sample, it is not 
expected to have the same ^-distribution as the underly- 
ing population. However, we can check that the observed 
sample is consistent with the assumed Galaxy model. Using 
the best-fit power law input <E>, we construct a smooth cu- 
mulative probability distribution function for z, from many 
"detected" CV samples. We then use a Kolmogorov-Smirnov 
(KS) test to compare this to the ^-distribution of the real 
sample of 20 systems; this is done many times, in order to 
sample the large errors in z of the real CVs. We find that 
the probability that the model and observed distributions 
are drawn from the same parent population is 0.56, for our 
model scale-height of 260 pc. Given that the observed sam- 
ple is small, this is not a very sensitive test — the sample 
is also consistent with an assumed scale height of 120 pc 
(in this case, the probability that the model and observed 
distributions are drawn from the same population is 0.14). 
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1O g(P0,R BS /P C " 3 ) 




-5 

lo g(Po,NEp/P C ~ 3 ) 

Figure 6. The probability distribution functions of log(po/pc) 
from the RBS (upper panel) and NEP data (lower panel). The 
vertical dashed lines mark 1-cr intervals. 



5.3 Consistency of the CV samples from the two 
surveys 

In iPretorius et~afl (|2007bJ) . we found p = 1.1± ;| x 
10" 5 pc" 3 from the NEP sample of only 4 CVs. Given the 
large errors in both this previous measurement and the one 
based on the combined survey presented here, they are not 
too different. However, if we calculate po based on the RBS 
alone, the result is almost an order of magnitude less than a 
measurement based only on the NEP survey. The log(po/pc) 
probability distribution functions for the two separate sam- 
ples are shown in Fig. [6] and are at first sight inconsistent. 

This would be easy to understand if the NEP survey 
were simply more powerful than the RBS, i.e. if it reached a 
fainter population (represented by the faint system EX Dra) 
than the RBS was capable of finding. However, this is not 
the case — the RBS is in fact more powerful than the NEP 
survey. Despite the much brighter flux limit of the RBS, it 
reached a larger volume (at all Lx) than the NEP survey, 
because of its wider angle (see Fig. [7]). It is then a fluke 
that the faintest CV in our combined sample was detected 
in the NEP survey, rather than in the RBS. Based on the 
detection of EX Dra in the NEP survey, we would expect to 
detect around 2 such systems in the RBS (since the volume 
of the RBS for CVs as bright as EX Dra is roughly twice 
that of the NEP survey); the non-detection of any CV at 
that Lx in the RBS is therefore unlucky, but has less than 
2-cr significance. 

Beyond this one dominant system, to determine 
whether the results of the two surveys can be reconciled, we 
need to consider the luminosities of all the observed systems. 
Fig. [8] shows the observed luminosity functions (constructed 
in the same way as for the combined survey in Section[T2]) of 




28 



30 

log(L x /erg s" 



32 



Figure 7. The logarithm of the survey volume of the RBS (fine 
curve), the NEP survey (dotted curve), and the combined survey 
(bold curve), as a function of Lx- This is basically the inverse 
of what is plotted in Fig. 1. The inset shows volume on a linear 
scale, for a small range in Lx; at the luminosity of EX Dra, 
log(Lx/erg s — 1 ) = 29.7, the RBS volume is roughly twice that of 
the NEP survey. 



-6 



-10 



RBS 
NEP 



29 



30 31 
log(L x /erg s" 1 ) 



32 



Figure 8. The observed and simulated <J> for the RBS (bold his- 
togram and large points) and NEP (fine histogram and smaller 
points) samples separately, where the simulated $ in both cases 
uses the best- fit input power law 4>, found for the combined survey 
(i.e., a = —0.80). Although the NEP and RBS observed luminos- 
ity functions are quite different, each is satisfactorily fit by the 
corresponding simulated "observed" for this assumed underly- 
ing luminosity function. 
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the RBS and NEP separately. Over-plotted are the output 
distributions, found as is Section 15.2.11 but again treating 
the two surveys separately. Both simulations use as input 
the best-estimate power law $ found in the previous sec- 
tion. Least-squares fits of the output from the models to the 
observed luminosity functions give reduced x 2 of 0.8 and 0.9 
for the RBS and NEP survey, respectively; both acceptable 
values. Therefore, despite the large difference in po when the 
two surveys are considered separately, they are consistent to 
within their uncertainties. 



6 DISCUSSION 

A basic assumption of the 1/V m ax method is that the de- 
tected objects are representative of the luminosity function 
of the true underlying CV population. Note that 'represen- 
tative' here does not imply that faint systems are as common 
in the observed sample as intrinsically, but only that some 
(or even one) are detected (see Section [3. 2. 3 [I . For an indi- 
cation of whether our sample is likely to be representative 
of the underlying CV population, and thus gives reliable p 
and $ estimates, it is important to know how faint CVs can 
be in X-rays. 

Theory predicts that the vast ma j ority of CVs shoul d 
be intrinsically f aint (e .g. IKolbl Il993l; iHowell et all Il997h : 
iPretorius efatl i|2007al ) and iPretorius fc Kniggei i|2008bl ) 
find that, although an as yet undetected faint CV popu- 
lation cannot dominate the overall population to the extent 
predicted, observed CV samples are nevertheless strongly 
biased against faint systems. The faintest secular average 
Lx expected of CVs can be estimate d from the gravitational 
radiatio n-driven M . We find that the lPatterson fc Ravmondl 
i|l985al ) relation between Lx and M predicts that the ma - 
jority of CVs in the theoretical population of IKolbl (Il993l ) 
should have time-averaged X-ray luminosities of a few times 
10 29 ergs -1 and higher. In intrinsically faint CVs, however, 
the rate of transfer of material onto the white dwarf surface 
(which determines Lx), is not the the same as the secular 
M, since these systems are dwarf novae. It is possible that, 
in the faintest CVs, hardly any material reaches the white 
dwarf surface during quiescence, so that they may spend 
most of their time at very faint Lx (perhaps with X-ray 
emission from the donor star being brighter than from the 
accretion flow). 

The faintest short-period CVs in our sample have 
luminosities of a few times 10 30 ergs -1 . Several intrinsi- 
cally faint, short-perio d systems are now k nown to have 
L x < 10 29 ergs -1 s fsee iBvckling et alJlioiol . as well as Pe- 



ter Wheatley, public communication^) . It is not known how 
intrinsically common such systems are, but they may well 
dominate the population. The standard theory of CV evo- 
lution pred icts that ~70% of CVs are period bouncers (e.g. 
IKolbl 1 19931 ) Using the observed mass-radius relationship of 
CV donors, iKnigge et al.l (|20 1 if ) predict an even larger frac- 
tion of period bouncers. Observations are also now indicat- 
ing a large population of intrinsically faint CVs (probably 
both normal short-period CVs and period bouncers). First, 



8 In a presentation at the conference "Wild 
Stars in the Old West", published on-line at 



http: / / www.noao.edu/meetings/wildstars2 / talks/ Wednesday /wheatley jnidsnatppt. 



iGansicke et alj (|2009h find a large number of intrinsically 
faint CVs at the shortest orbital periods. Furthermore, sev- 
eral period bouncers and good cand i date period bouncers 
are now known (e.g. iLittlefair et alj 20061; Littlefair et al.l 
|200S| ; iPattersonl 120111 ). and iPattersonl l|201ll) argue that 
these systems may be common enough to make up most 
of the intrinsic population. 

Two (related) properties of the sample used here show 
that it probably does not fairly represent the underlying 
population: it contains no faint short-period CVs, and it 
probably contains no period bouncers (see Section [2]). The 
lack of period bouncers alone likely means that it has 
missed at least half the intrinsic population. Furthermore, 
it is disconcerting that the faintest member of our sample, 
and therefore the system that dominates our p measure- 
ment, is a long-per iod CV, EX Drjfl Population synthe- 
sis models predict that at most a few percent of all CVs 
are a bove the period ga p (|Kolbl [l993 finds less than 1%, 
while IKnigge et al.|[201ll predict 3%). Although we find that 
long-period systems account for slightly more than 50% of 
our total space density, the data do not rule out these the- 
oretic al predictions. For example, using the IKnigge et al.1 
Hon]) fraction of long-period systems, and assuming that 
we have not significantly under-estimated the space density 
of long-period CVs, the space density of short-period CVs 
is ~ 2 x 10 -6 pc" 3 (97/3) ~ 6 x 10" 5 pc -3 . Using the upper 
limit on p from Section l5.ll we find that a short-period CV 
population of this size could escape detection in the two sur- 
veys, provided that the systems have Lx < 8 x 10 28 ergs -1 
(for the simple case of a hypothetical single-Lx population 
of faint, undetected CVs). 

Clearly, if CVs are arbitrarily faint in X-rays, the data 
allow for an arbitrarily large population to escape detec- 
tion. However, we can also choose to place a restriction 
in terms of what we might consider reasonable X-ray lu- 
minosities for active CVs. For example, if we integrate 
the best-fit power law luminosity function over the range 
28.7 < log(Lx/erg s _1 ) < 29.7 (this is a luminosity range 
where CVs are known to exist, but where we detect none), 
we find po = 1-2 x 10 -5 pc -3 for systems at those luminosi- 
ties, a factor of almost 3 larger than our po estimate from 
the detected CVs. This again indicates that it is reasonable 
to think that our po estimate is low by a factor of more than 
2. 

We find that a power law X-ray luminosity function ex- 
tending beyond log(_Ljf/erg s _1 ) ~ 32 is inconsistent with 
th e observed sample. Th i s is in agreement with figure 7 
of iPatterson fc Ravmondl (|l985aT ) , which shows that non- 
magnetic CVs should not be found at X-ray luminosities 
above ~ 10 32 ergs -1 . 

The main difficulty in constructing a luminosity func- 
tion from the RBS and NEP sample is that the distance 
estimates are very imprecise (with errors > 50% in some 
cases, implying that Lx is only very poorly constrained; see 



9 This is also discussed in IPretorius et alj|2007bl. where it was 
a more severe problem, since the sample was smaller. EX Dra 
accounts for more than half of the total po, and is almost en- 
tirely responsible for the faint end of the luminosity function. 
This system has Lx — 5 x 10 29 ergs -1 , which is unusually low 
for a long-P or b CV, and is probably caused by very high orbital 
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Section [2}. Since parallax measurements are at the moment 
only available for a fairly small number of relatively bright 
CVs, this problem will persist until results from a mission 
such as Gaia are available. 

We finally note that we have not considered the impact 
of variability on the X-ray luminosity function. The RASS 
data were taken over a period of about 6 months, during 
which the ecliptic poles were observed many ti mes, whereas 
lower ecliptic latitudes were covered only once l|Voges et al.l 
1999). Frequently outbursting DNe in the NEP area were 
therefore probably observed both in quiescence and in out- 
burst, while most DNe in the RBS sample were likely ob- 
served only in quiescence. 



7 CONCLUSIONS 

To summarize, we have constructed a complete, purely X- 
ray flux-limited sample of 20 non-magnetic CVs, and have 
used it to place constraints on the space density and X-ray 
luminosity function of CVs. Our main conclusions are listed 
below. 

(i) With the assumption that the combined non-magnetic 
CV sample from the RBS and NEP surveys is representative 
of the intrinsic population, we find po = 4+t x 10" 6 pc" 3 . 

(ii) It is likely that this po estimate excludes a large 
population of faint CVs (consisting of both normal short- 
period systems and period bouncers), and that it is low 
by at least a factor of ~2, as a result. In other words, the 
data are consistent with more than half of all CVs having 
28.7 < log(Lx/ergs _1 ) < 29.7, and escaping detection. 

(iii) To reach po — 2 x 1CP 4 pc~ 3 (at the high end of the 
predicted range), we require that the majority of CVs have 
L x < 4 x 10 28 ergs _1 . 

(iv) The precision with which p and $ can be measured 
is mainly limited by poorly constrained distances to CVs. 

(v) We find it impossible to correct for bias in a measure- 
ment of the X-ray luminosity function, without assuming a 
functional form for $. 

(vi) If the X-ray luminosity function of non-magnetic CVs 
is a truncated power law, $ = const x L\, the power law 
index that best reproduces our data is a = —0.8. 
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